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ELECTROMAGNETIC  THEORY  OF  REVERBERATION  CHAMBERS 


David  A.  Hill 

Radio  Frequency  Technology  Division 

National  Institute  of  Standards  and  Technology 

Boulder,  Colorado  80303 


This  report  presents  the  electromagnetic  theory  of  reverberation  chambers  as 
applied  to  electromagnetic  compatibility  measurements.  For  radiated  immunity 
measurements,  mode  theory  and  the  plane-wave  integral  representation  are  useful 
in  determining  electric  and  magnetic  field  characteristics  and  chamber  quality 
factor  (0.  The  plane-wave  integral  representation  also  provides  a  convenient 
method  for  describing  the  response  of  a  reference  antenna  or  a  test  object  in  a 
stirred  field.  For  radiated  emissions,  energy  conservation  or  electromagnetic 
reciprocity  can  be  used  to  describe  the  measurement  of  total  radiated  power. 
Comparisons  of  theory  and  measurements  are  presented  for  chamber  Q, 
probability  density  functions  of  fields  and  received  power,  and  radiated  immunity 
and  emissions  of  a  simple  test  object.  Recommendations  are  made  for 
improvements  in  statistical  electromagnetic  theory  as  applied  to  mechanical 
stirring  and  the  statistics  of  immunity  and  emissions  measurements. 

Key  words:  electromagnetic  compatibility;  emissions;  immunity;  mode-stirred 
chamber;  mode  theory;  ray  theory;  reverberation  chamber;  statistical 
electromagnetics. 


1.  INTRODUCTION 

The  use  of  reverberation  chambers  (also  called  mode-stirred  chambers)  for 
electromagnetic  compatibility  (EMC)  measurements  was  first  proposed  in  1968  [1].  It  took 
some  time  for  reverberation  chamber  measurements  to  gain  acceptance,  but  their  use  has 
increased  rapidly  in  the  past  decade.  The  practical  aspects  of  radiated-immunity  testing  in 
reverberation  chambers  were  thoroughly  covered  in  1986  in  NBS  Technical  Note  1092  [2],  and 
that  publication  has  since  served  as  the  primary  reference  handbook. 

Because  significant  advances  in  reverberation  chamber  theory  [3-5]  and  measurements 
[6-8]  have  been  made  since  the  publication  of  Technical  Note  1092  in  1986,  there  is  a  need  to 
update  the  information  in  that  document.  The  amount  of  relevant  new  material  has  become  so 
extensive  that  NIST  has  decided  to  split  it  into  three  reports.  The  purpose  of  this  report  is  to 
present  the  electromagnetic  theory  of  reverberation  chambers,  both  in  mathematical  detail  and  in 
physical  explanation.  The  remaining  two  publications  cover  statistical  analysis  [9]  and  practical 
applications  [10]. 

Reverberation  chambers  are  electrically  large,  high-(?  cavities  that  obtain  statistically 
uniform  fields  by  either  mechanical  stirring  [2]  or  frequency  stirring  [11,12].  This  report  will 


cover  only  mechanical  stirring  because  it  is  the  more  widely  used  method  and  frequency  stirring 
is  not  applicable  to  emissions  measurements.  Because  mechanically  stirred  reverberation 
chambers  are  large,  complex  cavities,  a  number  of  theoretical  methods  are  needed  to  provide  an 
electromagnetic  characterization. 

The  organization  of  this  report  is  as  follows.  Section  2  covers  cavity  mode  theory  which 
provides  useful  information  on  mode  fields,  mode  excitation,  and  mode  density,  all  of  which 
impact  the  effectiveness  of  mechanical  stirring.  Section  3  covers  a  plane-wave,  integral 
representation  which  allows  the  plane-wave  coefficients  to  be  random  variables  [13]  for 
describing  well-stirred  fields  in  reverberation  chambers.  This  statistical  field  theory  is  useful  in 
deriving  chamber  Q,  electric  and  magnetic  field  properties,  and  responses  of  antennas  and  test 
objects.  Reverberation  chambers  are  reciprocal  devices,  and  Section  4  covers  power 
conservation  and  reciprocity  theory  applied  to  radiated  emissions  measurements.  Conclusions 
and  recommendations  for  further  study  are  covered  in  Section  5.  Appendix  A  shows  how  ray 
theory  is  related  multiple  image  theory  which  can  provide  guidelines  for  mechanical  stirrer 
geometries.  Appendix  B  gives  the  theory  for  short  electric  dipoles  which  are  idealized  models 
for  electric-field  probes.  Appendix  C  gives  the  closely  related  theory  for  small  loop  antennas 
which  are  idealized  models  for  magnetic-field  probes.  Appendix  D  relates  chamber  Q  to 
chamber  decay  time. 


2.  MODE  THEORY 

The  theory  of  electromagnetic  cavities  with  separable  geometries  is  well  developed  in  the 
literature  [14-16],  but  the  application  of  the  theory  to  reverberation  chambers  is  not  straight- 
forward. Reverberation  chambers  typically  have  a  large,  complex  stirrer  (or  multiple  stirrers) 
and  operate  as  multimode  (rather  than  single-mode)  cavities.  The  stirrer  makes  it  difficult  to 
compute  the  chamber  fields,  even  for  a  simple  source  in  a  rectangular  cavity.  Also,  the  stirrer  is 
typically  rotated  through  a  large  number  of  positions  so  that  a  reverberation  chamber  is  really  an 
ensemble  of  a  large  number  of  cavities  with  different  shapes.  A  simpler  case  where  the  stirrer  is 
replaced  by  lateral  movement  of  one  wall  of  a  rectangular  cavity  has  been  analyzed  [17],  but  this 
geometry  is  not  used  in  practice  because  wall  movement  is  difficult  and  the  resultant  mode 
stirring  is  too  limited.  Stirrers  have  also  been  called  tuners  [2]  to  distinguish  between  stepped 
motion  and  continuous  motion,  but  the  term  stirrer  will  be  used  throughout  this  report. 

Despite  the  complexity  of  real  reverberation  chambers,  the  theory  of  empty  cavity  modes 
[18]  is  useful  in  predicting  some  important  chamber  properties,  such  as  mode  density  and  quality 
factor  (Q).  Also,  the  empty  cavity  modes  form  a  complete  set  for  expanding  the  fields  of  the 
same  cavity  when  it  contains  source  currents  (on  the  transmitting  antenna)  and  secondary 
scattering  currents  (on  the  stirrer)  [19].  This  section  will  concentrate  on  rectangular  cavities 
because  that  is  the  geometry  of  most  reverberation  chambers. 

2. 1  Empty  Cavity  Modes 

Consider  a  closed  cavity  of  arbitrary  shape  with  perfectly  conducting  walls  as  shown  in 
figure  1 .  The  interior  medium  (usually  free  space)  has  permittivity  s  and  permeability  fi.  The 
electric  field  Ep  of  thepth  eigenmode  satisfies  the  vector  Helmholtz  equation  where  the 
wavenumber  k  of  the  medium  takes  on  discrete  eigenvalues  kp  [14]: 


{v2+k;)Ep=0.  (1) 

The  electric  field  also  satisfies  the  divergence  equation 

V«£,=0.  (2) 

On  the  cavity  wall  the  tangential  electric  field  is  zero: 

hxEp=0,  (3) 

where  A  denotes  a  unit  vector. 

Without  loss  of  generality,  each  of  the  electric  eigenvectors  Ep  is  chosen  to  be  real.  The 
corresponding  magnetic  eigenvector  Hp  is  determined  from  Maxwell's  curl  equation 

//    =^Vx£d)  (4) 

where  tj  =  ^jju/s  .  Thus  Hp  is  purely  imaginary.  Physically,  the  electric  and  magnetic  fields  are 
standing  waves  in  phase  quadrature  oscillating  harmonically  with  an  angular  frequency  cop. 

:-t£-  (5) 

vV 

The  total  time-averaged  energy  Up  in  a  cavity  mode  is  equal  to  the  sum  of  the  time-averaged 
electric  energy  Upe  and  magnetic  energy  Upm  which  are  equal  [14]: 


Up=Upe+Upm     and     Upe=Upm.  (6) 

The  time-averaged  electric  and  magnetic  energies  are  defined  as  volume  integrals  [20]: 

u~=\\\l^iv  and  ^=f§WdF-  (7) 

Z      V  Z      V 

(The  coefficients  of  the  integrals  in  eq  (7)  involve  factors  of  XA  rather  than  V*  because  we  follow 
Harrington's  notation  [20]  of  using  RMS  values  rather  than  peak  values  for  phasors.) 

Numerical  methods  are  required  to  solve  eqs  (1)  through  (3)  for  the  eigenvalues  and 
eigenvectors  of  cavities  of  general  shape.  For  separable  geometries,  such  as  rectangular,  circular 
cylindrical,  or  spherical,  the  eigenfunctions  and  eigenvalues  are  well  known  [21].  Spherical  and 
circular  cylindrical  cavities  are  not  of  much  interest  for  reverberation  chambers  because  the 
curved  surfaces  can  produce  caustics  (focusing  of  rays)  that  make  it  difficult  to  obtain  spatial 
field  uniformity.  Complex  cavity  shapes  have  been  studied  numerically  in  the  physics 
community  to  study  chaotic  mode  fields  [22],  and  the  relevance  of  this  work  to  reverberation 


chambers  remains  an  open  question.  In  this  report,  we  will  concentrate  on  the  rectangular  cavity 
because  it  is  the  primary  geometry  in  use,  it  has  been  shown  to  produce  good  field  uniformity 
with  stirring,  and  it  is  the  best  understood  cavity  shape. 

The  geometry  for  an  empty  rectangular  cavity  is  shown  in  figure  2.  There  is  no  preferred 
axis,  and  all  cavity  modes  can  be  written  as  either  transverse  electric  (TE)  or  transverse  magnetic 
(TM)  to  any  of  the  three  coordinate  axes.  To  follow  the  usual  convention  in  the  literature  [18- 
20],  we  choose  to  construct  modes  that  are  either  TE  or  TM  to  the  z  direction.  The  electric  fields 
of  the  TE  and  TM  modes  can  be  written  in  the  following  forms  [19]: 

Ef(r)  =  -xkyrp{r)  +  yhxyp{r)  (8) 

and 

2?,™  =  -xkjctf  (r)  -  ykyhzpp  (r)  +  £(**  +  k]  )fp  (r) ,  (9) 

where 

Vpir)  =  -jJ=cos(kxx)sm(kyy)sm(k2z), 


VP{r)  =  -^sm(kxx)cos(kyy)sm(kzz),  (10) 


£ 

#(r)  =  ~HH  sin(***)  sin(Ar^)  cos(k2z), 
yjabc 


and 


In    .     _  mn        _nn 

kx  —      ,  k   —  — — ,  kz  —       , 
a  b  c 

(11) 

(2     for  7  =  0,  m  =  0,orn  =  0 

Sp~\'Jz     forl,m,n*0. 

The  modal  subscript/?  is  shorthand  for  the  nonnegative  integers  /,  m,  and  n.  Only  one  of  these 
three  integers  can  equal  zero  for  any  one  mode.  The  time  dependence,  exp(-/'ft^/),  is  suppressed. 
At  this  point,  the  normalization  of  the  modes  is  arbitrary,  but  the  above  normalization  is  useful  in 
deriving  the  Green's  function  for  the  source  problem  [19]. 

The  expressions  for  the  electric  field  of  the  TE  and  TM  modes  in  eqs  (8)  and  (9)  are 
consistent  with  other  results  in  the  literature  [18,20,21],  and  the  corresponding  magnetic-field 
expressions  can  be  obtained  by  substituting  eq  (8)  or  (9)  into  eq  (4).  Equation  (8)  shows  that  the 


z  component  of  the  electric  field  is  zero  for  TE  modes.  Equations  (4)  and  (9)  show  that  the  z 
component  of  the  magnetic  field  is  zero  for  TM  modes. 

2.2  Fields  of  Current  Sources 

In  this  section  we  consider  the  fields  excited  by  a  source  current  distribution  J(r' )  in  a 
rectangular  cavity  as  shown  in  figure  3.  The  current  and  the  induced  fields  have  a  time 
dependence  exp(-/&tf)  which  is  suppressed  in  the  rest  of  this  report.  The  induced  electric  field 
E{r)  can  be  written  as  an  integral  over  the  source  volume: 

E(r)  =  -icojuj\j    G(r,r>).J(r')dV,  (12) 

v, 

where  G  is  the  dyadic  Green's  function  which  satisfies  the  differential  equation  [19] 

V xV y.G{r,r')-k2G{r,r>)  =  18{r  -/•'),  (13) 


/is  the  identity  dyad,  5 is  the  three-dimensional  delta  function,  and  k  =  coy/jus  is  the 
wavenumber  of  the  medium.  The  dyadic  Green's  function  also  satisfies  the  boundary  condition, 

hxG(r,r')  =  0,  (14) 

for  either  r  or  r'  on  the  cavity  wall. 

The  dyadic  Green's  function  can  be  constructed  from  the  source-free  cavity  modes  in  eqs 
(8)  and  (9)  and  a  set  of  functions  that  have  nonzero  divergence  and  zero  curl  [19].  This  set  of 
functions  is  necessary  in  order  satisfy  Maxwell's  equations  in  the  source  region.  The  result  for 
G  [19]  can  be  written  in  terms  of  the  <f>  functions  which  were  defined  in  eq  (10): 

/=0  m=0  *=0  K    \Kp        K    ) 

~  %KkyP,  (r)fp  (/" )  -  xzKKfi  (r)fp  (/"  ) 

-yxkykJyp(rWp(r')  +  yy(k2  -k2yWp{rWP{r<)  (15) 

"  yzkykz^(rWp(r' )  -  zxkzkj;(r)<f>;(r' ) 

~ WW,W;P )  +  m2  - k)Wp{rWp{r<)}. 

This  result  is  in  agreement  with  the  dyadic  Green's  function  in  a  more  recent  publication  [17]. 

Several  comments  regarding  the  form  of  G  in  eq  (15)  are  in  order.  Each  of  the  nine 
dyadic  terms  in  the  triple  summation  satisfies  the  wall  boundary  condition  in  eq  (14),  so  G 
automatically  satisfies  eq  (14).  The  convergence  of  the  triple  sum  in  eq  (15)  becomes  very  slow 
as  the  observation  point  r  approaches  the  source  point  r' .  In  fact  the  sums  do  not  converge  when 
r  =  r'  because  G  is  singular.  The  details  of  the  singularity  of  G  and  the  evaluation  of  the  electric 
field  from  eq  (12)  have  been  thoroughly  studied  [23],  and  the  convergence  of  the  sums  can  be 


improved  by  converting  eq  (15)  to  a  sum  of  rays  and  modes  [19].  The  relationship  between  rays 
and  modes  will  be  discussed  further  in  Appendix  A. 

Because  of  the  factor  (k2p  -  k2)~] ,  eq  (15)  also  has  singularities  at  all  of  the  resonant 

frequencies  determined  by  k  =  kp  or  co  =  cop .  These  singularities  occur  only  for  lossless  cavities 
where  the  quality  factor  Q  is  infinite.  For  finite,  but  large,  values  of  cavity  Q,  an  approximate 

correction  [14]  to  eq  (15)  can  be  obtained  by  replacing  kp  by  kp{\ )  or  cop  by  cop{\ ) . 

This  replacement  keeps  the  modal  terms  finite  for  real  values  of  k  or  co  and  yields  a  half-power 
modal  bandwidth  of  copIQ.  If  Q  is  mode  dependent,  then  the  modal  Qp  should  be  used  in  the 
above  replacements,  and  the  half-power  modal  bandwidth  is  C0plQp. 

The  polarization  properties  of  fields  produced  by  currents  within  the  cavity  (on  either 
source  antenna  or  the  stirrer)  can  be  determined  from  eq  (15)  since  it  is  written  in  nine  separate 
dyads.  For  example,  consider  the  fields  produced  by  z-directed  currents,  J(r' )  =  zJz(r'),  as 
would  be  supported  by  a  z-directed  linear  antenna.  Then  eq  (12)  would  become 


E{r)  =  -icon\\\G{rS)*Uz(fW\  (16) 


The  vector  properties  of  the  electric  field  are  determined  by  the  following  dot  product  that 
extracts  three  of  the  nine  terms  in  G: 


co        co       co 


G(r,r').z=zZ2:  k,J  k>  H*,M;(r»;(o 

;=0m=0n=0       K     \Kp        K     )  (17) 

-  ykykzVP  {rWp (/"  )  +  z(k2  -  k)  Wp  (r)fp  (/•' )}. 

The  dot  product  in  eq  (17)  includes  all  three  vector  components;  so  the  electric  field  in  eq  (16) 
will  also  include  all  three  vector  components.  Even  though  z-directed  currents  will  excite  all 
three  components  of  the  electric  field,  they  do  not  provide  a  thorough  mode  excitation  of  the 
cavity.  The  way  to  see  this  is  to  derive  the  magnetic  field  from  Maxwell's  curl  equation: 

H  =  —VxE.  (18) 

ikrj 

If  eqs  (16)  and  (17)  are  substituted  into  eq  (18),  it  is  found  that  Hz  =  0.  This  is  because  z-directed 
currents  do  excite  TM  modes  of  the  form  given  by  eq  (9),  but  do  not  excite  TE  modes  of  the 
form  given  by  eq  (8).  Either  x-directed  or  ^-directed  currents  will  excite  TE  modes. 

Normally  an  effective  asymmetrical  stirrer  will  support  all  three  components  of  electric 
current  and  will  excite  both  TE  and  TM  modes.  The  task  of  computing  currents  on  electrically 
large,  complex  stirrers  is  extremely  difficult  (particularly  for  a  large  number  or  orientations);  so 
the  currents  needed  to  evaluate  the  electric  field  in  eq  (12)  are  generally  unknown.  Some  two- 
dimensional  stirrer  calculations  have  been  done  [24],  but  numerical  methods  are  only  beginning 
to  be  applied  to  realistic  three-dimensional  cavities  with  large  stirrers.  However,  the  rectangular- 


cavity  Green's  function  can  be  used  to  analyze  frequency  stirring  in  rectangular  cavities  [25] 
where  no  mechanical  stirrer  is  present. 

2.3  Mode  Number  and  Density 

The  source-free  modes  of  an  empty  rectangular  cavity  have  been  thoroughly  studied  by 
Liu,  Chang,  and  Ma  [18],  and  their  results  for  mode  number  and  density  are  summarized  in  this 
section.  For  the  rectangular  cavity  shown  in  figure  2,  the  eigenvalues  (wavenumbers)  kp  are 
determined  from 


2         /  \2         /  n.  2 

'  nn  » 


( ln\      (  mn^ 


\a  J 


Kn=®Lv£=  —    +  -r    +  —    ,  (19) 


\  c  J 


where  /,  m,  and  n  are  nonnegative  integers.  The  shorthand  subscript  p  will  be  replaced  by  the 
full  designation  Imn  in  this  section  to  make  the  mode  identification  more  explicit. 

In  mode  counting,  the  mode  degeneracies  need  to  be  considered.  When  none  of  the 
indicies  Imn  is  zero,  there  are  two  types  of  modes,  TE/W„  and  TM/m„,  for  each  eigenvalue.  So  the 
mode  degeneracy  is  2.  When  /  =  0,  only  a  TEom«  mode  exists  for  each  eigenvalue.  When  m  =  0, 
only  a  TEio„  mode  exists  for  each  eigenvalue.  When  n  =  0,  only  a  TM/wo  mode  exists  for  each 
eigenvalue.  The  forms  of  the  electric  fields  of  the  TE  and  TM  modes  are  given  in  eqs  (8)  and 
(9).  For  mode  counting,  we  designate  N{k)  as  the  number  of  modes  with  eigenvalues  k\m„  less 
than  or  equal  to  k.  This  total  number  of  modes  is  the  sum  of  the  number  of  modes  of  each  of  the 
above  types: 

N(k)  =N](k)  +  N2(k)  +  N3(k)  +  N4  (k)  +  N5(k) 

(20) 
=  2^I(*)  +  ^3(ifc)  +  ^4(ifc)  +  ^5(*), 

where  Ni  represents  TM/OT„  modes,  N2  (=  N\)  represents  TE/W„  modes,  A3  represents  TM/wo 
modes,  JV4  represents  TE0mn  modes,  and  N5  represents  TE/on  modes.  The  notation  and  properties 
of  the  five  mode  types  are  summarized  in  table  1. 

Equation  (20)  can  be  evaluated  numerically  by  mode  counting,  but  simple  analytical 
approximations  are  also  available.  The  actual  mode  number  is  a  discontinuous  function  of  £,  but 
the  best  smooth  approximation  Ns  is  [18] 

Ns{k)^e-a+b+ck+1-.  (21) 

The  first  term  on  the  right  side  of  eq  (21)  is  Weyl's  classical  approximation  Nw,  which  is  valid 
for  cavities  of  general  shape  and  can  be  written  in  terms  of  the  volume  V: 

Vk3 
NfV(k)  =  V^.  (22) 


The  extra  terms  in  eq  (21)  are  specific  to  the  rectangular  shape.  The  mode  numbers  in  eqs  (21) 
and  (22)  can  also  be  written  as  functions  of  frequency/: 


— abc^-ta  +  b  +  c)  —  +  — 
3  v3  v      2 


K(f)  =  —  abc^-(a  +  b  +  c)^-  +  -  (23) 


and 


tf,C/)=-2£4.  <24> 

3     v 

where  v  is  the  speed  of  light  in  the  medium  (usually  free  space).  Equations  (21)  through  (24)  are 
asymptotic  high-frequency  approximations  that  are  valid  when  the  cavity  dimensions  are 
somewhat  greater  than  a  half  wavelength. 

Numerical  results  for  N  (by  computer  counting),  Ns,  and  Nw  are  shown  in  figure  4  for  the 
NIST  reverberation  chamber  {a  =  2.74  m,  b  =  3.05  m,  and  c  =  4.57  m).  The  extra  terms  in  Ns 
improve  the  agreement  obtained  with  Weyl's  formula.  The  smooth  mode  density  Ds(f)  is  also 
shown  in  figure  4.  It  is  obtained  by  differentiating  eq  (23): 


<W,(/)     „      ,    f     a  +  b  +  c 
d/  v3 


DXf)  =  ^LL  =  *7tabc±T .  (25) 


The  Weyl  approximation  again  equals  the  first  term: 

Dw(f)  =  ^P-  =  ^V^-.  (26) 

d/  v3 

The  mode  density  is  an  important  chamber  design  parameter  because  it  determines  how  many 
modes  are  present  in  a  small  bandwidth  about  a  given  frequency.  For  example,  figure  4  shows 
that  the  NIST  chamber  has  a  mode  density  somewhat  greater  than  1  mode  per  megahertz  at  a 
frequency  of  200  MHz.  Experience  has  shown  that  the  NIST  chamber  provides  adequate 
performance  at  frequencies  above  200  MHz,  but  not  below  200  MHz  where  the  mode  density  is 
too  low  to  obtain  spatial  field  uniformity  [2]. 


3.  PLANE-WAVE,  INTEGRAL  REPRESENTATION 

Deterministic  mode  theory,  as  discussed  in  the  previous  section,  is  useful  in  determining 
empty-cavity  properties,  but  it  is  not  convenient  for  predicting  the  response  of  a  receiving 
antenna  or  test  object  in  a  mechanically  stirred  reverberation  chamber.  Since  many  stirrer 
positions  are  employed  in  reverberation  chamber  measurements,  some  type  of  statistical  method 
[4]  is  required  to  determine  the  statistics  of  antenna  or  test  object  response.  At  the  same  time,  it 
is  important  to  ensure  that  the  associated  electromagnetic  theory  is  consistent  with  Maxwell's 
equations. 


The  purpose  of  this  section  is  to  present  a  plane-wave,  integral  representation  for  the 
electric  and  magnetic  fields  that  satisfies  Maxwell's  equations  and  also  includes  the  statistical 
properties  expected  for  a  well-stirred  field  [13].  The  statistical  nature  of  the  fields  is  introduced 
through  the  plane-wave  coefficients  that  are  taken  to  be  random  variables  with  fairly  simple 
statistical  properties.  Because  the  theory  uses  only  propagating  plane  waves,  it  is  fairly  easy  to 
use  to  calculate  the  responses  of  test  objects  or  reference  antennas. 

Section  3.1  presents  the  basic  plane- wave  integral  theory  and  derivations  of  important 
field  properties.  Section  3.2  includes  derivations  of  the  responses  of  antennas  or  test  objects  to 
the  statistical  field.  Section  3.3  derives  probability  density  functions  for  the  fields  and  test  object 
responses.  Section  3.4  analyzes  the  loss  mechanisms  and  the  Q  of  a  reverberation  chamber. 

3 . 1  Electric  and  Magnetic  Field  Properties 

A  typical  geometry  for  an  immunity  measurement  in  a  reverberation  chamber  is  shown  in 
figure  5.  A  transmitting  antenna  radiates  cw  fields,  and  the  mechanical  stirrer  (or  multiple 
stirrers)  is  rotated  to  generate  a  statistically  uniform  field.  The  test  volume  can  occupy  a  fairly 
large  portion  of  the  chamber  volume. 

The  electric  field  E  at  location  r  in  a  source-free,  finite  volume  can  be  represented  as  an 
integral  of  plane  waves  over  all  real  angles  [26] 


E(r)  =  jjF(Q)exp(/A  •  r)dO,  (27) 


4m 


where  the  solid  angle  Q  is  shorthand  for  the  elevation  and  azimuth  angles,  a  and  /?,  and 

dQ.  =  sin  a  da  d/?  .  The  exp(-/a#)  time  dependence  is  suppressed.  The  vector  wavenumber  k  is 

k  =  -£(jcsin<zcos/?  +  j>  sin  a  sin /?  +  £  cos  a) .  (28) 

The  angular  spectrum  F(Q)  can  be  written 

F(Q)  =  aFa(Q)  +  ^(Q),  (29) 

where  a  and  f3  are  unit  vectors  that  are  orthogonal  to  each  other  and  to  k.  Both  Fa  and  Fp  are 
complex  and  can  be  written  in  terms  of  their  real  and  imaginary  parts: 

Fa(Q)  =  Far(Cl)  +  iFfir(a)    and    F/}(Q)  =  F^ty  +  iF^Q).  (30) 

The  geometry  for  a  plane-wave  component  is  shown  in  figure  6. 

The  electric  field  in  eq  (27)  satisfies  Maxwell's  equations  because  each  plane-wave 
component  satisfies  Maxwell's  equations.  For  a  spherical  volume,  the  representation  in  eq  (27) 
can  be  shown  to  be  complete  because  it  is  equivalent  to  the  rigorous  spherical-wave  expansion 
[27].  For  a  non-spherical  volume,  the  plane-wave  expansion  can  be  analytically  continued 
outward  from  a  spherical  volume,  but  the  general  conditions  under  which  the  analytic 


continuation  is  valid  have  yet  to  be  established.  In  this  report,  we  assume  that  the  volume  is 
selected  so  that  eq  (27)  is  valid. 

For  a  statistical  field  as  generated  in  a  reverberation  chamber,  the  angular  spectrum  F(Q) 
is  taken  to  be  a  random  variable  (which  depends  on  stirrer  position).  For  derivation  of  many  of 
the  important  field  quantities,  the  probability  density  function  of  the  angular  spectrum  is  not 
required,  and  it  is  sufficient  to  specify  certain  means  and  variances.  In  a  typical  reverberation 
chamber  measurement,  the  statistical  ensemble  is  generated  by  rotating  the  stirrer  (or  stirrers).  In 
this  report,  <  >  represents  an  ensemble  average  over  stirrer  position.  The  starting  point  for  the 
statistical  analysis  is  to  select  statistical  properties  for  the  angular  spectrum  that  are 
representative  of  a  well-stirred  field  which  would  be  obtained  in  an  electrically  large,  multimode 
chamber  with  a  large,  effective  stirrer  [2].  Appropriate  statistical  assumptions  for  such  a  field 
are  as  follows: 

<Fa(Q)>  =  <Ffi(Q)>  =  0,  (31) 

<FQr(n])^(n2)>  =  <i>(Q1)F(a2)>  = 

<Far(Ql)Ffir(Cl2)>=<Far(a,)F^Q2)>=  (32) 

<Fai(Qx)Ffir(Cl2)>  =  <Fai(Cll)Ffi(n2)>=0, 

and 

<Far(Q])Far(Q2)>^<Fai(Q,)Fai(Q2)>  = 

(33) 
<^r(Q1)^r(n2)>  =  </v,(^1)iv,(Q2)>=C^(Q1-Q2), 

where  5  is  the  Dirac  delta  function  and  Ce  is  a  constant  with  units  of  (V/m)  . 

The  mathematical  reasons  for  the  assumptions,  eqs  (31)  through  (33),  will  become  clear 
when  the  field  properties  are  derived,  but  the  physical  justifications  are  as  follows.   Since  the 
angular  spectrum  is  a  result  of  many  rays  or  bounces  with  random  phases,  the  mean  value  ought 
to  be  zero,  as  indicated  in  eq  (31).  Since  multipath  scattering  changes  the  phase  and  rotates  the 
polarization  many  times,  angular  spectrum  components  with  orthogonal  polarizations  or 
quadrature  phase  ought  to  be  uncorrelated,  as  indicated  in  eq  (32).  Since  angular  spectrum 
components  arriving  from  different  directions  have  taken  very  different  multiple  scattering  paths, 
they  ought  to  be  uncorrelated  as  indicated  by  the  delta  function  on  the  right  side  of  eq  (33).  The 
coefficient  Ce  of  the  delta  function  is  proportional  to  the  square  of  the  electric  field  strength  as 
will  be  shown  later.  The  following  useful  relationships  can  be  derived  from  eqs  (32)  and  (33): 

<Fa(Q,)/v(Q2)>=0,  (34) 

and 
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<Fa(Ql)F:(Q2)>  =  <Ffi(Q,)F;(Q2)>  =  2CES{Q]  -Q2),  (35) 

where  *  denotes  complex  conjugate. 

A  number  of  field  properties  can  be  derived  from  eqs  (27)  and  (31)  through  (35). 
Consider  first  the  mean  value  of  the  electric  field  <E>,  which  can  be  derived  from  eqs  (27)  and 
(31): 

<  E{r)  >  =  {{  <  F(Q)  >  exp(/*  •  r)  dO.  =  0 .  (36) 


An 


Thus  the  mean  value  of  the  electric  field  is  zero  because  the  mean  value  of  the  angular  spectrum 
is  zero.  This  result  is  expected  for  a  well-stirred  field  which  is  the  sum  of  a  large  number  of 
multipath  rays  with  random  phases. 

The  square  of  the  absolute  value  of  the  electric  field  is  important  because  it  is 
proportional  to  the  electric  energy  density  [5].  From  eq  (27),  the  square  of  the  absolute  value  of 
the  electric  field  can  be  written  as  a  double  integral: 

\E(r)\2=j\jJF(n]).F\Q2)exp[i(k,-k2)T]dQ,dn2.  (37) 

An  An 

The  mean  value  of  eq  (37)  can  be  derived  by  applying  eqs  (34)  and  (35)  to  the  integrand: 

<|£(r)|2>=  4CEjjj\d(Q,  -Q2)exp[z(^  -*2).r]dQ,dQ2 .  (38) 

An  An 

One  integration  in  eq  (38)  can  be  evaluated  by  using  the  sampling  property  of  the  delta  function, 
and  the  second  integration  is  easily  evaluated  to  obtain  the  final  result: 

<\E(r)\2>  =  4CEjjdn2  =  \6nCE=E20.  (39) 

An 

Thus  the  mean-square  value  of  the  electric  field  is  E2  and  is  independent  of  position.  This  is  the 

field  uniformity  property  of  an  ideal  reverberation  chamber;  it  applies  to  the  ensemble  average  of 
the  squared  electric  field  and  has  been  verified  experimentally  with  an  array  of  electric-field 
probes  [2], [8].  For  convenience  throughout  the  remainder  of  this  report,  Ce  is  defined  in  terms 
of  the  mean-square  value  of  the  electric  field  as  indicated  in  eq  (39). 

By  a  similar  derivation,  the  mean-square  values  of  the  rectangular  components  of  the 
electric  field  can  be  derived: 

<\Ex\2>  =  <\Ey\2>  =  <\Ez\2>=^-.  (40) 

This  is  the  isotropy  property  of  an  ideal  reverberation  chamber,  and  it  has  been  verified  with 
three-axis,  electric-field  probes  [2], [8]. 
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The  magnetic  field  H  can  be  derived  by  applying  Maxwell's  curl  equation  to  eq  (27): 

H(r)  =  —  VxF(r)  =  -ff*xF(Q)exp(/*»r)dQ.  (41) 

Applying  eq  (3 1)  to  eq  (41)  shows  that  the  mean  value  of  the  magnetic  field  is  zero: 

<  H(r)  >  =  -  f f  k  x  <  F(Q)  >  exp(/A  •  r)  dQ  =  0 .  (42) 

The  square  of  the  magnitude  of  the  magnetic  field  can  be  written 

\H(r)\2=  Ajjjjft  xF(Q,)].[a2  xF*(Q2)]exp[/(^  -*2).r]dQ,  dQ2 .  (43) 


4n  4n 


The  derivation  of  the  mean-square  value  follows  closely  that  of  the  electric  field,  and  the  result  is 

<\H(r)\2>  =  ^.  (44) 

n 

Thus  the  mean-square  magnetic  field  also  exhibits  spatial  uniformity,  and  the  value  is  related  to 
the  mean-square  electric  field  by  the  square  of  the  free-space  impedance: 

<\H(r])f>=<E{rl)f>,  (45) 

rj 

where  r\  and  r-i  are  arbitrary  locations.  This  free-space  relationship  has  been  demonstrated 
experimentally  by  using  electric  and  magnetic  field  probes  [2]. 
The  energy  density  W  can  be  written  [20] 

W(r)  =  ±[s\E(r)\2+ju\H(r)\2].  (46) 

The  mean  value  can  be  obtained  from  eqs  (39),  (44),  and  (46): 

<  W{r)  >  =  ^[e<\  E(r)  f  >  +M  <\  H(r)  \2 >]  =  aE2 .  (47) 

Thus  the  average  value  of  the  energy  density  is  also  independent  of  position. 
The  power  density  or  Poynting  vector  S  can  be  written  [20] 

S(r)  =  E(r)xH\r).  (48) 
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From  eqs  (27),  (42),  and  (48),  the  mean  of  the  power  density  can  be  written 

<  S(r)  >  =  -  j]  JJ  <  F(Q,  )x[k2xF*  (Q2 )]  >  exp[/(*,  -  k2 )  •  r]  dQ,  dQ2 .  (49) 

I  An  Alt 

The  expectation  in  the  integrand  can  be  evaluated  from  vector  identies  and  eqs  (34)  and  (35): 

F2 

< F(Q])x[k2xF(Q2)]>  =  k2 -^£(Q,  -Q2).  (50) 

An 

The  right  side  of  eq  (49)  can  now  be  evaluated  from  eq  (50)  and  the  sampling  property  of  the 
delta  function: 

E2    rr~ 
< S(r)  >=  -=2-  \\k2dQ2  =  0.  (51) 

A  physical  interpretation  of  eq  (51)  is  that  each  plane  wave  carries  equal  power  in  a  different 
direction  so  that  the  vector  integration  over  An  steradians  is  zero.  This  result  is  important 
because  it  shows  that  power  density  is  not  the  proper  quantity  for  characterizing  field  strength  in 
reverberation  chambers.  The  mean  value  of  energy  density  as  given  by  eq  (47)  is  an  appropriate 
positive  scalar  quantity  that  could  be  used.  Another  possibility  is  to  define  a  positive  scalar 
quantity  S  that  has  units  of  power  density  and  is  proportional  to  the  mean  energy  density: 

E2 
S  =  v<W>==2.,  (52) 

1 

where  v  =  1  /  -^/ue  .  For  lack  of  a  better  term,  S  will  be  called  scalar  power  density  in  the  rest  of 

this  report.  This  quantity  could  be  used  to  compare  with  uniform-field,  plane-wave  testing  where 
power  density,  rather  than  field  strength,  is  sometimes  specified. 

To  this  point,  field  properties  at  a  point  have  been  considered.  Real  antennas  and  test 
objects  have  significant  spatial  extent,  and  the  spatial  correlation  function  [28]  of  the  fields  is 
important  in  understanding  responses  of  extended  objects  in  reverberation  chambers  [29].  The 
spatial  correlation  function  pir\,r2)  of  the  electric  field  can  be  defined  as 

<E(r) •  E*(n)  >  ,„x 

4<\E(rx)\2x\E{r2)\2> 

where  r\  and  r2  are  two  arbitrary  locations.  The  numerator  of  eq  (53)  is  the  mutual  coherence 
function  which  has  been  used  to  describe  wave  propagation  in  random  media  [30].  The 
denominator  of  eq  (53)  can  be  evaluated  from  eq  (39),  and  the  numerator  can  be  evaluated  from 
eqs  (27),  (34),  and  (35)  so  that  the  final  result  is  [28] 
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P(r„r,)= 'fa<*"'-'i  P.  (54) 

k\rx-r2  | 

The  identical  correlation  function  has  been  derived  from  cavity  mode  theory  [4]  and  radiative 
transfer  theory  [31]  and  has  been  checked  experimentally  [29].  The  same  correlation  function 
can  be  derived  for  the  magnetic  field,  and  it  also  applies  to  acoustic  reverberation  chambers  [32]. 
A  correlation  length  lc  can  be  defined  as  the  separation  corresponding  to  the  first  zero  in  eq  (54): 

klc=7r     or    lc=7r/k  =  X/2,  (55) 

where  X  is  the  wavelength  in  the  medium. 

An  angular  correlation  function  p(s,,  s2)  can  be  defined  as 

,-    -x  <ErAr)E.(r)>  ,r^ 

p(*.,0=  ,  Y  ,   »  <56) 

V<|£il(r)|2x|£,2(r)|2> 
where  the  two  electric  field  components  are  defined  as 

Esl(r)  =  s^E(r)    and     Es2(r)  =  s2  .  E(r)  (57) 

and  s,  and  s2  are  unit  vectors  separated  by  an  angle  /as  shown  in  figure  7.  From  eq  (40),  the 
denominator  of  eq  (56)  is  found  to  equal  El  13 .  The  numerator  of  eq  (56)  is  evaluated  from  eqs 
(27),  (34),  and  (35),  and  the  result  for  the  angular  correlation  function  is 

p(s1,s2)  =  s^s2=cosy.  (58) 

This  result  is  independent  of  position  r.  The  same  angular  correlation  function  can  be  derived 
for  magnetic  field  components.  Equation  (58)  shows  that  the  three  rectangular  components  of 
the  electric  field  are  uncorrelated,  and  this  is  consistent  with  the  theory  of  Kostas  and  Boverie 
[33]. 

3.2  Antenna  or  Test-Object  Response 

Consider  now  a  receiving  antenna  or  a  test  object  placed  in  the  test  volume.  The  simplest 
case  of  a  lossless,  impedance-matched  antenna  will  be  considered  first.  The  received  signal  can 
be  written  as  an  integral  over  incidence  angle  by  analogy  with  Kerns's  plane-wave,  scattering- 
matrix  theory  [34].  The  received  signal  could  be  a  current,  a  voltage,  or  a  waveguide  mode 
coefficient,  but  the  general  formulation  remains  the  same.  Consider  the  received  signal  to  be  a 
current  /  induced  in  a  matched  load.  For  an  antenna  located  at  the  origin,  the  current  can  be 
written  as  a  dot  product  of  the  angular  spectrum  with  a  receiving  function  Sr(Q)  integrated  over 
angle: 
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7  =  JJ5r(Q)«F(Q)dQ,  (59) 


An 


where  the  receiving  function  can  be  written  in  terms  of  two  components, 

sr(n)  =  asra(n)+flsrfi(ci).  (60) 

In  general,  Sra  and  Srp  are  complex,  so  the  antenna  can  have  arbitrary  polarization,  such  as  linear 
or  circular.  For  example,  a  z-directed  linear  antenna  with  linear  polarization  would  have 
Srp(Q)  =  0  •  A  circularly  polarized  antenna  would  have  Sr/}(Q)  -  ±iSra(Q)  for  right-  or  left- 
hand  circular  polarization. 

The  mean  value  of  the  current  I  can  be  shown  to  be  zero  from  eqs  (31)  and  (59): 

< I  >=  j"jX(Q)«  <  F(Q)  >  dQ  =  0.  (61) 

All 

The  absolute  value  of  the  square  of  the  current  is  important  because  it  is  proportional  to  received 
power  Pr: 

Pr  =|/|2  Rr  =/?rJJJJ[5r(ni).F(QI)][5;(n2).F*(Q2)]dQ1dn2,  (62) 

Ax  An 

where  the  radiation  resistance  Rr  of  the  antenna  is  also  equal  to  the  real  part  of  the  matched  load 
impedance.  The  mean  value  of  the  received  power  can  be  determined  from  eqs  (34),  (35),  and 
(62): 

<Pr>-<\I\2>Rr=^-^j\[\Sra(Cl2)\2  +  \Srfi(Q2)\2]dn2.  (63) 


2  Ax 


Art 


The  physical  interpretation  of  eq  (63)  is  that  the  ensemble  average  of  received  power  is  equal  to 
an  average  incidence  angle  (Q2)  and  polarization  {a  and  /?  components). 

The  integrand  of  eq  (63)  can  be  related  to  the  effective  area  of  an  isotropic  antenna  A2/4;r 
and  the  antenna  directivity  D(Q.2)  by  [35] 


^      £ 


rjRr[\Sra(^2)r  +IW)    ]  =  - D(Q2).  (64) 

An 

Substitution  of  eq  (64)  into  eq  (63)  yields 

<Pr>  =  L%j?L±UD{Cl2)dCl2.  (65) 

2  j]  4n  4n  J[ 
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The  integral  in  eq  (65)  is  known  because  the  average  (over  D2)  of  D  is  1.  Thus  the  final  result 
for  the  average  received  power  is 

<P>=±f^^_.  (66) 

2   77   4/r 

The  physical  interpretation  of  eq  (66)  is  that  the  average  received  power  is  the  product  of  the 
scalar  power  density  E\lr]  times  the  effective  area  X2 I  An  of  an  isotropic  antenna  times  a 

polarization  mismatch  factor  of  one  half  [36].  This  result  is  independent  of  the  antenna 
directivity  and  is  consistent  with  the  reverberation  chamber  analysis  [37]  of  Corona  et  al.  Some 
of  the  earlier  data  indicated  that  eq  (66)  was  in  better  agreement  with  measurements  if  the  one 
half  polarization  mismatch  factor  was  omitted  [2].  However,  more  recent  comparisons  of 
antenna  received  power  with  field-probe  data  [8]  and  with  a  well-characterized  test  object  [38] 
support  the  inclusion  of  the  factor  of  one  half.  Consequently,  the  polarization  mismatch  factor 
needs  to  be  included  to  be  in  agreement  with  theory  and  with  most  measured  data.  Traditionally, 
linearly  polarized  antennas  have  been  used  as  reference  antennas  in  reverberation  chambers,  but 
this  analysis  suggests  that  circularly  polarized  antennas  are  also  appropriate.  Experimental  data 
with  circularly  polarized  would  be  useful  for  confirming  this  theoretical  result.  The  special  cases 
of  an  electrically  short  dipole  (electric-field  probe)  and  an  electrically  small  loop  (magnetic-field 
probe)  are  discussed  in  Appendices  B  and  C. 

The  preceding  analysis  can  be  extended  to  the  case  of  a  real  antenna  with  loss  and 
impedance  mismatch  by  using  Tai's  theory  [36].  The  effective  area  Ae  can  be  generalized  to 

At(Cl)  =  ^-D(Cl)pmrfat  (67) 

Ak 

where/?  is  the  polarization  mismatch,  m  is  the  impedance  mismatch,  and  rja  is  the  antenna 
efficiency.  All  three  quantities,  p,  m,  and  rja,  are  real  and  can  vary  between  0  and  1 .  The 
average  of  Ae  over  incidence  angle  and  polarization  can  be  written  [37] 

A2 
<Ae>=-—mrja.  (68) 

The  average  received  power  is 

<Pr>=^o_<Ae>,  (69) 

where  E\  I  rj  can  again  be  interpreted  as  the  average  scalar  power  density. 

Test  objects  can  be  thought  of  as  lossy,  impedance-mismatched  antennas,  so  eq  (69)  also 
applies  to  test  objects  as  long  as  terminals  with  linear  loads  can  be  identified.  This  theory  has 
been  used  to  predict  the  responses  of  an  apertured  coaxial  line  [39],  an  apertured  rectangular  box 
[5],  and  a  microstrip  transmission  line  [38,40]  as  compared  to  a  reference  antenna  in  a 
reverberation  chamber.  Good  agreement  with  measurements  has  been  obtained  in  each  case. 


16 


The  microstrip  example  is  a  good  illustration  of  the  use  of  the  above  theory.  The 
response  of  terminated  microstrip  transmission  line  was  computed  using  the  above  theory  and 
measured  in  the  NIST  reverberation  chamber  [40]  using  the  setup  shown  in  figure  8.  A 
comparison  of  theory  and  measurements  is  shown  in  figure  9  for  frequencies  from  200  to  2000 
MHz.  The  plotted  quantity  is  the  ratio  of  the  average  power  received  by  the  reference  antenna  to 
the  average  power  received  by  the  microstrip  line  in  decibels.  (This  ratio  is  sometimes  called 
shielding  effectiveness  in  decibels.)  The  theoretical  ratio  is  201og]0[(A2  /8;r)/  <  Ae  >],  where 

/l2/8;ris  the  theoretical  average  effective  area  of  the  reference  antenna  and  <A<>  is  the  average 
effective  area  of  the  microstrip  transmission  line.  The  measurements  were  performed  on  three 
different  physical  models,  and  the  "bottom  feed"  microstrip  line  best  fits  the  theoretical  model. 
Hence  it  has  the  best  agreement  with  theory.  Even  that  measured  curve  has  a  small  negative  bias 
which  is  probably  due  to  impedance  mismatch  in  the  reference  antenna  which  was  not  taken  into 
account.  The  actual  reference  antenna  was  a  log  periodic  dipole  array  below  1000  MHz  and  a 
broadband  ridged  horn  above  1000  MHz. 

3.3  Probability  Density  Functions 

The  statistical  assumptions  for  the  angular  spectrum  in  eqs  (31)  through  (33)  have  been 
used  to  derive  a  number  of  useful  ensemble  averages  in  sections  3.1  and  3.2.  These  results  have 
not  required  a  knowledge  of  the  particular  form  of  the  probability  density  functions.  However, 
such  knowledge  would  be  very  useful  for  analysis  of  measured  data  which  is  always  based  on 
some  limited  number  of  samples  (stirrer  positions). 

The  starting  point  for  deriving  electric-field  probability  density  functions  is  to  write  the 
rectangular  components  in  terms  of  their  real  and  imaginary  parts: 

Ex=Exr+iExi,    E^E^+iE^,    Ea=Ew+iE<.  (70) 

(The  dependence  on  r  will  be  omitted  where  convenient  because  all  of  the  results  in  this  section 
are  independent  of  r.)  The  mean  value  of  all  the  real  and  imaginary  part  in  eq  (70)  is  zero,  as 
shown  in  eq  (36): 

<Exr  >  =  <£„.  >  =  <Eyr>  =  <Eyi  >  =  <Ezr  >=<Ezi  >=0.  (71) 

The  variance  of  the  real  and  imaginary  parts  can  be  shown  to  equal  half  the  result  for  the 
complex  components  in  eq  (40): 

<E2xr>  =  <El>  =  <E2yr  >  =  <£*  >=<E*,>  =  <E*>=&  =  a\  (72) 

o 

The  mean  and  variance  of  the  real  and  imaginary  parts  in  eqs  (71)  and  (72)  are  all  the 
information  that  can  be  derived  from  the  initial  statistical  assumptions  in  eqs  (31)  and  (33). 

However,  the  maximum  entropy  method  [41], [42]  can  be  used  to  derive  the  probability 
density  functions  from  eqs  (71)  and  (72).  For  example,  consider  Exr.  The  maximum  entropy 
method  selects  the  probability  density  function  J[Exr)  to  maximize  the  entropy  (uncertainty)  given 
by  the  integral 
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-  ]f(EJWf(Exr)]dExr  (73) 

-00 

subject  to  the  constraints  in  eqs  (71)  and  (73)  and  the  usual  probability  constraint, 

]f(Exr)dExr=\.  (74) 

-00 

The  maximization  of  eq  (73)  subject  to  the  constraints  in  eqs  (71),  (72),  and  (74)  is  performed  by 
the  method  of  Lagrange  multipliers  [42],  and  the  result  for/  is  the  normal  distribution: 


*J27T<J 


2a2 


(75) 


where  a  is  defined  in  eq  (72).  The  same  probability  density  function  also  applies  to  the  other 
real  and  imaginary  parts  of  the  electric  field  components.  Even  though  there  are  other 
probability  density  functions  that  would  satisfy  the  constraints  in  eqs  (71),  (72),  and  (74),  the 
normal  distribution  in  eq  (75)  maximizes  the  entropy  given  by  eq  (73)  and  is  the  least  biased. 
Any  other  probability  density  function  would  have  to  be  based  on  additional  information  that  is 
not  provided  by  this  analysis. 

Equations  (27),  (33),  and  (34)  can  be  used  to  show  that  the  real  and  imaginary  parts  of  the 
electric-field  components  are  uncorrelated.  Only  the  derivation  for  <ExrExl>  will  be  shown,  but 
the  derivations  for  the  other  correlations  are  similar.  From  eqs  (27)-(3 1),  the  real  and  imaginary 
parts  of  Ex  can  be  written: 


Exr(r)  =  ^{[cosacospFar(Q)-sm^Ffir(Q)]cos(kT) 


4ff 


(76) 


-  [cos a  cos/?  Fa  (Q)  -  sin  /?  Fp%  (Q)]sin(A  •  r)}  dQ, 
Exi  (r)  =  jj  { [cos  a  cos  /?  Fa  (Q)  -  sin  /?  Fp>  (Q)]  cos(*  •  r) 

An 

+  [cos a  cos/?  Far(Q)  -  sin  /?7^r(Q)]sin(A  •  r)}dQ. 


(77) 


The  average  value  of  the  product  of  eqs  (76)  and  (77)  can  be  evaluated  by  using  eqs  (32)  and 
(33)  inside  the  double  integral  and  making  use  of  the  delta  function  to  evaluate  one  integration. 
Then  the  remaining  integrand  is  zero: 


E, 

< 


EAr)Exi(r)>=  -f- \\ [cos2  a2  cos2  /?2  +sin2  /?2][cos(/r2  •r)sin(A2  t) 

16;r  J/K  (78) 

-  cos(£2  •  r)sin(£2  •  r)]dQ  -  0. 


is 


Similar  evaluations  show  that  the  real  and  imaginary  parts  of  all  three  rectangular  components 
are  uncorrected.  Since  they  are  Gaussian,  they  are  also  independent  [43]. 

Since  the  real  and  imaginary  parts  of  the  rectangular  components  of  the  electric  field  have 
been  shown  to  be  normally  distributed  with  zero  mean  and  equal  variances  and  are  independent, 
the  probability  density  functions  of  various  electric  field  magnitudes  or  squared  magnitudes  are 
chi  or  chi-square  distributions  with  the  appropriate  number  of  degrees  of  freedom.  The 
magnitude  of  any  of  the  electric  field  components,  for  example  \EX\,  is  chi  distributed  with  two 
degrees  of  freedom  and  consequently  has  a  Rayleigh  distribution  [43]: 


f(\Ex\)=l-^exp 


2a' 


(79) 


Figure  10  shows  a  comparison  of  eq  (79)  with  measured  data  taken  at  1  GHz  in  the  NASA  A 
Chamber  [8].  The  chamber  has  two  stirrers,  and  the  total  number  of  samples  (stirrer  positions)  is 
225.  The  data  were  taken  with  a  small  electric-field  probe  that  was  calibrated  at  NIST  [8].  The 
agreement  is  about  as  good  as  can  be  expected  with  225  samples. 

The  squared  magnitude  of  any  of  the  electric  field  components,  for  example  |£x|2,  is  chi- 
square  distributed  with  two  degrees  of  freedom,  and  consequently  it  has  an  exponential 
distribution  [44]: 


f(\Exf)  = 


1 


2a' 


■exp 


2a' 


(80) 


The  probability  density  functions  in  eqs  (79)  and  (80)  agree  with  Kostas  and  Boverie  [33].  They 
suggest  that  the  exponential  distribution  in  eq  (80)  is  also  applicable  to  the  power  received  by  a 
small,  linearly  polarized  antenna,  but  it  will  shown  later  that  the  exponential  distribution  applies 
to  the  power  received  by  any  type  of  antenna. 

The  total  electric  field  magnitude  \E\  is  chi  distributed  with  6  degrees  of  freedom  and  has 
the  following  probability  density  function  [3 3], [43]: 


/(l^l) 


|£|5 
^6XP 


2a' 


(81) 


Figure  1 1  shows  a  comparison  of  eq  (81)  with  measured  data  taken  under  the  same  conditions  as 
in  figure  10.  In  this  case  a  three-axis,  electric-field  probe  was  used  to  take  the  data  [8].  Again 
the  agreement  is  about  as  good  as  can  be  expected  for  225  samples. 

The  squared  magnitude  of  the  total  electric  field  is  chi-square  distributed  with  6  degrees 
of  freedom  and  has  the  following  probability  density  function  [43]: 


/(|£|2) 


E\A 
—  exP 
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2a' 


(82) 
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The  dual  probability  density  functions  for  the  magnetic  field  can  be  obtained  by  starting  with  the 
variance  of  the  real  or  imaginary  parts  of  one  of  the  magnetic  field  components,  for  example  Hxr: 


<H„> 


Gni 


2    ~UH  • 


(83) 


Then  the  dual  of  the  results  in  eqs  (79)  through  (82)  can  be  obtained  by  replacing  E  by  H  and  o 

by  crH. 

Similar  techniques  can  be  used  to  analyze  the  signal  received  by  an  antenna.  Equation 
(61)  shows  that  the  real  and  imaginary  parts  of  the  current,  Ir  and  /,,  have  zero  mean: 


</  >  =  </. .  >=0. 


(84) 


The  derivation  of  eq  (66)  can  be  modified  to  obtain  the  variance  of  the  real  and  imaginary  parts 
of  the  current: 


rl  rl  1        ^n     ^ 

<I:  >  =  <!:  >- 


4Rr   r]  4k 


=  o 


i  ■ 


(85) 


Since  only  the  mean  and  variance  of  the  the  real  and  imaginary  parts  of  the  current  are  known, 
the  maximum  entropy  method  is  again  applicable  for  determining  the  probability  density 
function.  The  result  is  again  the  normal  distribution  for  both  Ir  and  /,: 


Air) 


1 


IttOj 


exp 


2o' 


and     /(/,.)  = 


1 


■yjZxOj 


exp 


2o) 


(86) 


Equations  (32),  (33),  (59),  and  (60)  can  be  used  to  show  that  Ir  and  /,  are  uncorrected: 

<IrI,>=0.  (87) 

Equations  (84)  through  (87)  can  be  used  to  show  that  the  current  magnitude  \I\  is  chi  distributed 
with  two  degrees  of  freedom  (Rayleigh  distribution): 


o, 


ill 
2o) 


(88) 


The  current  magnitude  squared  |/|  is  chi-square  distributed  with  two  degrees  of  freedom 
(exponential  distribution): 


/(|/|2) 


1 
2o) 


exp 


2o) 


(89) 
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From  eq  (62),  the  received  power  Pr  is  proportional  to  l)  + 1] .  So  Pr  is  also  chi-square 
distributed  with  two  degrees  of  freedom  and  has  an  exponential  probability  density  function: 


/«)=s?ic«p 


2a)Rr 


(90) 


If  the  load  is  not  matched,  then  the  radiation  resistance  Rr  in  eq  (90)  is  replaced  by  the  load 
resistance  Ri  and  the  expression  for  o)  is  modified  from  that  given  in  eq  (85).  The  result  in  eq 
(90)  is  in  agreement  with  that  of  Kostas  and  Boverie  [33],  but  its  derivation  is  more  general.  It 
does  not  rely  on  the  received  power  being  proportional  to  the  field  at  a  point  (such  as  an  electric 
or  magnetic  field  probe)  and  is  valid  for  general  extended  antennas.  The  same  result  is 
applicable  for  general  test  objects  as  long  as  they  are  linear  and  have  identifiable  terminals  with 
linear  loads. 

The  result  in  eq  (90)  matches  experimental  data  in  reverberation  chambers  fairly  well  for 
a  variety  of  antennas  (dipoles,  horns,  and  log-periodic  dipole  arrays).  Figure  12  shows  a 
comparison  of  eq  (90)  with  data  taken  with  a  log-periodic  dipole  array  at  1GHz  in  the  NASA  A 
Chamber  [8].  The  agreement  is  somewhat  better  than  in  figures  10  and  1 1  because  the  number 
of  samples  was  increased  to  1024.  In  figures  10  through  12,  the  transmitting  antenna  was  a  log- 
periodic  dipole  array. 

3.4  Loss  Mechanisms  and  Cavity  Q 

In  eq  (39),  El  was  introduced  as  the  mean-square  value  of  the  electric  field  which  was 
shown  to  be  independent  of  position.  This  constant  can  be  related  to  the  power  Pt  transmitted 
and  the  chamber  Q  by  conservation  of  power  [5],  [44].  The  starting  equation  is  the  definition  of 
quality  factor  (0: 

Q  =  f,  (9.) 

where  U  is  the  energy  stored  in  the  cavity  and  Pd  is  the  power  dissipated.  Since  the  average 
energy  density  was  shown  to  be  independent  of  position  in  eq  (47),  the  stored  energy  can  be 
written  as  the  product  of  the  average  energy  density  and  the  chamber  volume  V: 

U=<W>V.  (92) 

For  steady-state  conditions,  conservation  of  power  requires  that  the  dissipated  power  Pd  equals 
the  transmitted  power  Pt.  Then  eqs  (47),  (91),  and  (92)  can  be  used  to  derive 

£o2="%  (93) 

coev 
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This  analysis  can  be  carried  further  to  relate  the  transmitted  power  to  the  power  received 
by  a  receiving  antenna  located  in  the  chamber.  If  eq  (93)  is  substituted  into  eq  (66),  the  average 
power  received  by  a  matched,  lossless  antenna  is  found  to  be 


\6tt2V 


<Pr>=T7^Pt  (94) 


Equations  (93)  and  (94)  show  the  importance  of  the  Q  enhancement  in  determining  the  field 
strength  or  the  received  power  in  the  chamber.  The  most  popular  method  of  measuring  Q  is 
based  on  the  solution  of  eq  (94)  for  Q. 

Equation  (95)  is  applicable  to  an  impedance-matched,  lossless  receiving  antenna,  but  dissipative 
or  mismatch  losses  can  be  accounted  by  modifying  the  effective  area  as  shown  in  eq  (68). 

The  calculation  of  chamber  Q  requires  that  all  losses  are  accounted  in  evaluating  Pd  in  eq 
(91).  A  theory  has  been  developed  for  including  the  following  four  types  of  loss  [5]: 

^=^1+^2+^3+^4.  (96) 

where  Pdi  is  the  power  dissipated  in  the  cavity  walls,  Pdi  is  the  power  absorbed  in  loading 
objects  within  the  cavity,  Pdi  is  the  power  lost  through  aperture  leakage,  and  Pd*  is  the  power 
dissipated  in  the  loads  of  receiving  antennas.  By  substituting  eq  (96)  into  eq  (91),  we  can  write 
the  following  expression  for  the  inverse  of  Q: 

Q-^Qf+Qf+Qj+Q?,  (97) 

where 

coU  o)U  o)U  (oU 

Qi=—,Q2  =  -^-, &  =  "T->  and  Q4  =  — .  (98) 

"d\  "dl  *dl  *d4 

The  four  loss  mechanisms  can  be  analyzed  as  follows.  Wall  loss  is  usually  dominant,  so  it  will 
be  covered  in  most  detail. 

For  highly  conducting  walls,  the  plane-wave  integral  representation  can  be  analytically 
continued  all  the  way  to  the  wall  surfaces,  and  the  reflected  fields  are  related  to  the  incident 
fields  via  plane- wave  reflection  coefficients  as  shown  in  figure  13.  Then  Pdi  in  eq  (98)  can  be 
evaluated  in  terms  of  the  wall  area  ,4  and  the  wall  reflection  coefficient.  The  resultant  expression 
to  be  evaluated  for  Q\  is  [45] 

IkV 
n  = tfL 1  (99) 

1      ,4<(l-|r|2)cos0>Q 
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where  <  >q  indicates  an  average  over  incidence  angle  and  polarization  and  T  is  a  reflection 
coefficient. 

The  reflection  coefficients  for  TE  (perpendicular)  polarization  TVe  and  TM  (parallel) 
polarization  Ttm  are  given  by  [26] 


_  /dwk  cos  0-ti<Jkl  -k2sin20 

te  —  r^ ^ ; — 


juwk cos 6  +  /Jyjkl  -k2  sin2  $ 
and 


(100) 


uk2  cosO  -  u^kJk2  -  k2  sin2  6 
*■**=— 2 W      ;  f         =  •   .    2      ,  (101) 

f&l  cos 9  +  fiJcJkl  -kzsm2e 


where  &w  =  co ^///w  (sw  +  /'aw  /  <y) ,  aw  is  the  wall  conductivity,  ew  is  the  wall  permittivity,  and  u.v 

is  the  wall  permeability.  To  account  equally  for  both  polarizations  in  eq  (99),  the  average 
quantity  can  be  written 


<(i-|r|2)cos0>Q=<[i-l(|rTC  |2  +|r™  |2)]cos0>Q 


(102) 


=  jil-klT^l2  +\YTU\2)]cosesm0d0. 


o        2 


For  |  kw  Ik  |»  1 ,  the  squares  of  the  reflection  coefficients  can  be  approximated 

ir„f»i-4^*R**-)co,g  (103) 


and 


4MJKc(kJ 


M\K\   cos# 
where  Re  indicates  real  part.  Substitution  of  eqs  (102)  through  (104)  into  eq  (99)  yields 

1\k    I2  V 
Q,  »        '    wl  ,  (105) 

4^rRe(*w) 

where  /^  =  /Av///. 
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Equation  (105)  does  not  require  that  the  walls  be  highly  conducting.  However,  if  the 
walls  are  highly  conducting  and  aw  /{coe^ )  »  1 ,  then  Q\  simplifies  to 

a^  (106) 

where  S  «  2/  yJo)^iwaw  .  This  is  the  usual  expression  for  reverberation  chamber  Q  for  the  case 

where  wall  losses  are  dominant.  A  related  derivation  has  employed  the  skin  depth 
approximation  from  the  start  followed  by  an  average  over  an  ensemble  of  plane  waves  [3].  For 
the  case  of  a  rectangular  cavity  where  the  modes  are  known,  it  has  been  derived  by  averaging  the 
modal  Q  values  for  modes  whose  resonant  frequencies  are  in  the  vicinity  of  the  excitation 
frequency  [18].  The  modal  derivation  for  rectangular  cavities  was  done  for  the  case  //r  =  1,  and  a 
correction  was  derived  for  rectangular  cavities  [18].  The  correction  is  important  only  at  low 
frequencies. 

If  the  cavity  contains  absorbers  (lossy  objects  distinct  from  the  walls),  the  absorption  loss 
Pai  can  be  written  in  terms  of  the  absorption  cross  section  <ra  [46]  which  is  generally  a  function 
of  incidence  angle  and  polarization: 

Pd2=S,<cTa>a.  (107) 

The  appropriate  average  is  over  4;rsteradians  and  both  (TE  and  TM)  polarizations  [44]: 

<^>Q  =  ^{j(^TC+^™)dn.  008) 

The  absorption  cross  section  in  eq  (108)  can  be  that  of  a  single  object  or  a  summation  for 
multiple  absorbers.  For  example,  for  M absorbers  <ca  >Q  is  replaced  by 

M 
<^>a=^,<^>Qt  (109) 

where  <  cra.  >  is  the  averaged  absorption  cross  section  of  the  rth  absorber.  From  eqs  (98)  and 
(107)  the  result  for  Q2  is  [44] 

2nV 

Q2  =  ■  (no) 

A<CTa  >Q 

The  formulation  for  leakage  loss  Pdi  is  similar  to  that  of  absorption  loss  because  apertures 
can  be  characterized  by  a  transmission  cross  section  <t/  [47].  However,  only  plane  waves  that 
propagate  toward  the  wall  aperture(s)  contribute  to  leakage  power.  So  the  expression  for  Qi  is 
modified  from  eq  (1 10)  by  a  factor  of  2  [44]: 
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Q3  =  •  (in) 

A  <<Jl  >n 


Also,  the  angular  average  is  over  2tu steradians  (0  <  6  <  n 1 2 ): 

1 


JJVrm+^mJ^-  (H2) 


<<T,    >r,  =  

An 


2x 


For  the  case  of  L  apertures,  <  a,  >n  in  eq  (1 12)  is  replaced  by  a  summation: 


L 

<0-/>Q=S<a/'   >Q»  (U3) 

i=l 


where  <  ati  >a  is  the  averaged  transmission  cross  section  of  the  /th  aperture.  For  electrically 

large  apertures,  <  a,  >0  is  independent  of  frequency  and  Q3  is  proportional  to  frequency.  For 

small  or  resonant  apertures,  the  frequency  dependence  of  Q3  is  more  complicated.  The  Q  of  a 
cavity  with  a  circular  aperture  has  been  studied  theoretically  and  experimentally  [44]. 

The  power  dissipated  in  the  load  of  a  receiving  antenna  is  covered  in  Section  3.2.  For  a 
lossless  receiving  antenna,  Pd4  can  be  written 

P«=^S,  (114) 

where  m  is  the  impedance  mismatch.  From  eqs  (98)  and  (114),  Q4  can  be  written 

a=^f.  (.15) 

mA 

If  there  are  multiple  receiving  antennas,  eqs  (1 14)  and  (115)  can  be  modified  accordingly.  For 
example,  if  there  are  N  identical  receiving  antennas,  P^  is  multiplied  by  N  and  Qa  is  divided  by 
N.  For  a  matched  load  (m  =1),  Q4  is  proportional  to  frequency  cubed.  This  means  that  Q4  is 
small  for  low  frequencies  and  is  the  dominant  contributor  to  the  total  Q  in  eq  (97).  The  effect  of 
antenna  loading  on  the  Q  of  reverberation  chambers  has  been  observed  experimentally  [11].  At 
high  frequencies,  Q4  becomes  large  and  contributes  little  to  the  total  Q. 

A  comparison  of  measured  and  calculated  Q  [44]  is  shown  in  figure  14  for  a  rectangular 
aluminum  cavity  of  dimensions  0.514  m  x  0.629  m  x  1.75  m.  The  Q  measurements  were 
performed  by  the  power  ratio  method  of  eq  (95)  and  decay  time  method  as  discussed  in 
Appendix  D.  Standard-gain,  Ku-band  horn  antennas  were  used  to  cover  the  frequency  range 
from  12  to  18  GHz.  The  measured  Q  values  fall  below  the  theoretical  Q,  but  the  agreement  is 
much  better  than  that  obtained  in  earlier  comparisons  [2].  The  decay-time  measurement  [7] 
generally  agrees  better  than  the  power-ratio  method  because  it  is  less  affected  by  antenna 
efficiency  and  impedance  mismatch. 
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A  second  comparison  of  theory  and  measurement  in  figure  15  shows  the  effect  of  loading 
the  cavity  with  three  spheres  of  radius  0.066  m  filled  with  salt  water  [44].  In  this  case  the 
absorption  loss  as  described  by  eq  (1 10)  decreases  the  Q  dramatically.  Broadband  ridged  horns 
were  used  for  the  Q  measurements,  and  the  agreement  with  theory  is  not  as  good.  However,  the 
decay-time  measurement  is  again  a  significant  improvement  over  the  power-ratio  measurment. 


4.  RADIATED  EMISSIONS 

Reverberation  chambers  have  been  used  primarily  for  radiated  immunity  measurements, 
and  as  a  result  a  great  deal  of  research  has  been  done  in  characterizing  chamber  fields.  However, 
reverberation  chambers  are  reciprocal  devices,  and  they  can  and  have  been  used  for  radiated 
emissions  measurements  [39].  The  quantity  measured  is  the  total  radiated  power,  and  the 
measurement  can  be  explained  by  either  power  conservation  [44]  or  reciprocity  [48]. 

4. 1  Radiated  Power 

If  the  equipment  under  test  (EUT)  radiates  (transmits)  power  Pieut,  eq  (94)  can  be  used 
to  determine  the  average  power  <PrEUT>  received  by  a  matched,  lossless  reference  antenna. 
Equation  (94)  is  based  on  conservation  of  power,  and  it  can  be  solved  for  Pieut- 


PtEVT   =  —f?T  <  PrEUT   >■  0  16) 


£Q 


In  theory  this  equation  could  be  used  directly  for  measurement  ofPtEUT-  However,  eq  (1 16) 
requires  that  the  chamber  volume  Fand  (loaded)  Q  be  known.  It  also  requires  that  the  receiving 
antenna  be  impedance  matched  and  lossless  or  that  the  received  power  be  corrected  for  antenna 
effects. 

A  better  way  to  determine  Pieut  is  to  perform  a  separate  reference  measurement  under  the 
same  chamber  conditions.  If  a  known  power  Ptref\s  transmitted  and  an  average  power  <Prref>  is 
received,  the  coefficient  on  the  right  side  of  eq  (1 16)  can  be  determined: 


\6k2V         P. 


tref 


(117) 


*Q  <  Prref   > 

Then  Pieut  can  be  determined  by  the  ratio 

p 

PtEVT   ~  7T <  PrEUT   >  ■  0  ]  8) 

<Prref> 

If  the  same  receiving  antenna  is  used  for  both  the  EUT  and  the  reference  measurement,  this 
method  has  the  additional  advantage  of  approximately  canceling  efficiency  and  impedance 
mismatch  effects  of  the  receiving  antenna. 
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This  was  done  in  the  measurement  of  radiated  power  from  a  microstrip  line  [40],  and  the 
agreement  between  theory  and  measurement  as  shown  in  figure  1 6  was  good.  The  actual 
quantity  plotted  was  the  following  power  ratio 

<P"ef:        Ptref  .  (119) 


<P       >     P 

Since  the  same  input  power  was  fed  to  the  reference  antenna  and  the  microstrip  line,  the  ratio  in 
figure  16  can  be  interpreted  as  either  a  shielding  effectiveness  or  the  reciprocal  of  the  radiation 
efficiency  of  the  microstrip  line. 

4.2  Reciprocity  Relationship  to  Radiated  Immunity 

Electromagnetic  reciprocity  has  many  mathematical  forms,  and  it  can  be  applied  to  fields, 
circuits,  or  a  mixture  of  the  two  [49].  Since  reciprocity  involves  interchanging  the  source  and 
receiver,  it  provides  a  method  for  relating  radiated  emissions  and  immunity.  Consider  an  EUT 
located  at  the  center  of  a  spherical  volume  as  shown  in  figure  17.  In  an  immunity  measurement, 
the  EUT  is  illuminated  by  incident  electric  and  magnetic  fields,  Et  and  //,,  due  to  sources  located 
outside  the  spherical  surface  Sr.  In  an  emissions  measurement,  the  EUT  radiates  (transmits) 
electric  and  magnetic  fields,  Et  and  Ht. 

A  typical  EUT  is  very  complex,  and  de  Hoop  and  Quak  [48]  have  developed  a  multi-port 
reciprocity  formulation  to  relate  emissions  and  immunity.  Here  we  consider  the  simpler  special 
case  of  a  single  port  within  the  EUT  as  shown  in  figure  18.  In  an  immunity  measurement  the 
incident  fields  induce  an  open-circuit  voltage  Vu  and  Z,  is  the  impedance  of  the  Thevenin 
equivalent  circuit.  An  arbitrary  load  impedance  Z/  is  connected  across  the  terminals.  In  an 
emissions  measurement,  Vt  is  zero  and  a  current  /,  flows  in  the  loop.  The  radiated  fields  are 
proportional  to  /,  and  can  be  normalized  as  follows: 

Et(r)  =  Ite„(r)    and     Ht(r)  =  I,hH(r),  (120) 

where  e„  and  h„  are  the  electric  and  magnetic  fields  that  are  radiated  when  /,  =  1  A.  If  reciprocity 
is  applied  at  the  circuit  terminals  and  the  spherical  surface,  the  following  expression  is  obtained 
for  V,  [48],[50]: 

Vi i  =-JJr*k(r)x^|.(r)-£I.(r)xAll(r)]d5r.  (121) 

Up  to  this  point,  eq  (121)  is  fairly  general  because  there  are  no  restrictions  on  the  sphere 
radius  r  or  the  incident  fields.  If  the  surface  integral  in  eq  (121)  is  performed  in  the  far  field  of 
the  EUT  (kr  »  1) ,  the  normalized  EUT  fields  can  be  written  in  the  following  forms 
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(  \        ra  ^\  exp(/Ar) 


h„(r)  =  rxet(8,t) 


(122) 
exp(/'Ar) 
rjr 


where  et{6,(f>)»r  =  0  and  # and  ^  are  standard  spherical  coordinates.  To  apply  eq  ( 1 2 1 )  to 

reverberation  chamber  measurements,  the  incident  electric  and  magnetic  fields  are  replaced  by 
plane-wave  integral  representations  from  eqs  (27)  and  (41).  Then  eq  (121)  can  be  rewritten  as 

1 A*  (123) 

-[ffF(Q)exp(/^«r)dQ]x[Ifxef(<9,^)]}dS,r. 

ii  n 

To  evaluate  the  surface  integration,  it  is  written  explicitly  in  terms  of  spherical 
coordinates: 

Inn 

jj{     }dSr  =  jj{     }r2sm0d0d</>.  (124) 

S,  0  0 

The  exponential  factor  exp(/A*r)  in  eq  (123)  is  a  rapidly  oscillating  function  of  6  and  <j>  except  at 
the  stationary  point  r  =  -k  .  A  stationary-phase  [51]  evaluation  of  eq  (123)  yields 

yi=?£L(U.  {et(a,p)  x  [k  x  F(Q)] 

kjl  H  (125) 

+  F(Q)x[kxe,(a,/])]}dCl. 

Because  the  reciprocity  integral  in  eq  (121)  is  independent  of  the  surface  over  which  it  is 
evaluated,  the  result  in  eq  (125)  is  an  exact,  rather  than  an  asymptotic,  result.  (This  is  consistent 
with  the  observation  that  eq  (125)  is  independent  of  r.)  Vector  identities  can  be  used  to  reduce 
eq(125)to 

yi=^LHel(aifi)*F(Cl)6n.  (126) 

kTi    An 

This  is  as  far  as  the  expression  for  V\  can  be  simplified.  It  shows  that  the  open-circuit 
voltage  induced  when  the  EUT  is  illuminated  in  an  immunity  test  is  proportional  to  a  weighted 
integral  of  the  transmitted  far  field  et  when  the  EUT  is  transmitting.  Equation  (126)  is  similar  to 
the  earlier  receiving  response  in  eq  (59)  except  that  the  receiving  function  in  eq  (59)  was  not 
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derived  in  terms  of  the  transmission  properties  of  the  antenna.  Another  interpretation  of  eq  (126) 
is  that  the  transmitting  and  receiving  patterns  of  an  antenna  or  an  EUT  are  the  same. 

The  statistical  properties  of  the  plane-wave  spectrum  F(Q)  were  discussed  in  Section  3.1, 
and  they  can  be  used  to  derive  the  statistical  properties  of  V\.  For  example,  eqs  (31)  and  (126) 
can  be  used  to  show  that  the  average  value  of  F,  is  zero: 

<  V1 >  =  —  \\ et{a,p)*  <  F(Q)  >dQ  =  0.  (127) 

kri  a 

The  mean  square  value  of  Vt  is  the  most  useful  quantity  because  it  is  proportional  to  the  received 
power  in  an  emissions  measurement.  The  squared  magnitude  \Vj\2  can  be  written 


VA2  = 


rA7u\ 


j  JJJjK(a1,A).F(Q1)]K,(a2,/?2).F,(n2)]dnidQ2.  (128) 


An  An 


The  average  value  <\Vi  |2>  can  be  determined  by  applying  the  properties  of  F  in  eqs  (34)  and 
(35)  to  eq  (128): 
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<\Vl\2>=^\\\el(a],^)\2dQ].  (129) 

k    *1       An 

From  eq  (122),  the  integral  on  the  right  side  of  eq  (129)  can  be  related  to  the  total  radiated  power 
Prad\  when  a  current  of  1  A  is  flowing  in  the  circuit  in  figure  14: 

P««  =~  jj \e„ (<*,fi) I2  dO  =  - \j \e, (a,  fi)  |2  dQ.  (130) 

"    An  'I  An 

Substitution  of  eq  (130)  into  eq  (129)  yields 


<W,t>  =  ^p-P^.  (131) 

Equation  (131)  shows  that  the  total  radiated  power  in  an  emissions  measurement  is 
proportional  to  the  mean-square,  induced  voltage  in  an  immunity  measurement.  For  an  arbitrary 
current  /  in  the  transmitting  (emissions)  case,  the  radiated  power  Prad  is  given  by 

Prad=\I\2Rrad,  (132) 

where  Rrad  is  the  radiation  resistance  part  of  the  transmitting  impedance  Z,  in  figure  18.  For 

/  =  1A ,  we  have  Pradx  -  Rrad .  If  we  substitute  for  Prad\  and  k  (=  2^/A),  eq  (13 1)  can  be  rewritten 
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|2^     /sad       \  12 


<i^r>/(4/^)_A' 


(133) 


The  numerator  of  the  left  side  of  eq  (133)  is  the  received  power  for  the  case  of  a  matched  load 
(ZL  =  Z*)  with  no  dissipative  loss  in  the  circuit  (Re(Z,)  =  Rrad)  in  figure  18,  and  the 
denominator  is  the  scalar  power  density.  This  ratio  is  the  average  effective  area,  and  it  is  equal 
to  A2/8;ras  shown  previously  in  eq  (66). 

If  the  circuit  in  figure  18  has  loss  (Re(Z,)  =  Rrad  +  Rhss)  but  is  still  impedance  matched 

(ZL  =  Z*) ,  eq  (133)  can  be  manipulated  to  the  following  form 


{<|  Vt  |2>  l[A{Rrad  +  Rloss)]}/{E20  In)  _       Rrad 


?i*x  Kt+Kss 


(134) 


In  eq  (134),  the  numerator  is  the  average  received  power  divided  by  the  scalar  power  density 
which  equals  the  average  effective  area.  The  denominator  X  /8;ris  the  maximum  effective  area 
for  any  receiving  antenna  in  a  well-stirred  field.  Kraus  [52]  has  termed  this  ratio  "effectiveness 
ratio,  a,"  for  the  simpler  case  where  the  incident  is  a  plane  wave  which  can  be  polarization 
matched  by  the  receiving  antenna  to  yield  a  maximum  effective  area  of  A2/4;r.  The  right  side  of 
eq  (134)  is  the  radiation  efficiency  r\a  for  the  emissions  case.  So  we  can  rewrite  eq  (134): 

ai  (immunity)  =  rja  (emissions).  (135) 

The  theoretical  and  experimental  results  in  figures  9  and  16  provide  a  verification  of  eq  (135)  for 
the  specific  case  of  a  microstrip  transmission  line.  Typically  in  the  EMC  community  the  left  side 
of  eq  (135)  is  called  shielding  effectiveness  and  is  given  in  decibels.  If  there  is  impedance 
mismatch,  both  sides  of  eq  (134)  can  be  multiplied  by  the  same  mismatch  factor  to  provide  a 
comparison  with  ideal  receivers  or  transmitters. 


5.  CONCLUSIONS  AND  RECOMMENDATIONS 

As  a  result  of  the  complexity  of  reverberation  chambers,  a  number  of  theoretical  methods 
are  needed  to  provide  an  electromagnetic  characterization.  Mode  theory  is  usually  the  method  of 
choice  for  cavity  analysis,  and  it  has  been  useful  for  determining  mode  density  and  modal  and 
composite  Q  for  rectangular  cavities  [18].  However,  it  is  not  convenient  for  quantitative  analysis 
of  real  chambers  with  mechanical  stirring  because  the  cavity  modes  are  hard  to  determine  and 
they  change  with  stirrer  rotation.  Ray  theory  can  provide  some  useful  insight,  but  it  is  also 
difficult  to  use  for  real,  high-Q  chambers  because  of  the  geometrical  complexity  and  the  slow 
convergence  of  ray  sums. 

A  plane-wave,  integral  representation  with  a  statistically  uniform  angular  spectrum  is 
useful  for  describing  well-stirred  fields.  It  has  the  advantage  of  being  easy  to  use  for  describing 
antenna  and  test-object  responses  because  it  uses  only  propagating  plane  waves.  It  is  most 
directly  useful  for  predicting  results  in  immunity  measurements  as  in  Section  3,  but  reciprocity 
can  be  used  to  predict  results  in  emissions  measurements  as  in  Section  4.  When  the  maximum 
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entropy  method  is  applied  to  the  results  of  the  integral  representation,  probability  density 
functions  can  be  derived  for  many  field  and  antenna  response  quantities.  These  statistical  results 
are  in  agreement  with  other  statistical  theories  that  are  based  on  the  central  limit  theorem  and 
with  experimental  results. 

A  number  of  theoretical  extensions  and  improvements  would  be  useful.  For  deterministic 
field  calculations,  modern  computational  tools  could  be  used  to  characterize  chamber  fields  and 
the  effectiveness  of  mechanical  stirring.  The  transmission  line  method  (TLM)  has  been  used  to 
study  a  simple  two-dimensional  model  of  mechanical  stirring  [24],  and  the  finite  element  method 
has  been  used  for  three-dimensional  analysis  of  chambers  [53].  The  use  of  the  rectangular  cavity 
Green's  function  in  an  efficient  form  [19]  has  been  suggested  for  method-of-moments  analysis  of 
mechanical  stirring.  Numerical  methods  are  most  likely  to  be  successful  at  lower  frequencies 
where  the  electrical  size  of  the  cavity  is  not  too  large,  and  this  is  an  important  frequency  range 
where  real  chambers  have  operational  problems  (such  as  ineffective  stirring,  poor  field 
uniformity,  and  antenna  impedance  mismatch).  A  long-standing  question  regarding  mechanical 
stirring  is  how  much  frequency  shift  of  a  given  mode  can  be  generated  by  stirrer  rotation  so  that 
a  meaningful  comparison  with  frequency  stirring  [1 1],[12]  can  be  made. 

Some  extensions  in  the  plane-wave,  integral  representation  would  also  be  useful.  The 
theory  is  most  applicable  to  ensemble  averages,  but  most  immunity  testing  uses  peak  values 
rather  than  average  values.  For  a  finite  number  of  samples,  the  probability  density  functions  in 
section  3.3  can  be  used  to  predict  peak  values,  and  this  has  worked  fairly  well  for  modest 
numbers  of  samples.  However,  these  peak  values  approach  infinity  as  the  number  of  samples 
approaches  infinity.  For  quantities  such  as  received  power,  this  is  not  physically  possible.  The 
plane-wave,  integral  representation  can  probably  still  be  used,  but  the  tails  of  the  distribution 
need  some  modification  to  satisfy  energy  conservation.  Another  idealization  that  requires  some 
further  attention  is  the  delta  function  for  angular  correlation  in  eq  (33).  For  real,  imperfect 
stirring,  the  delta  function  probably  needs  to  be  replaced  by  a  peaked  function  with  nonzero 
width.  A  final  point  that  requires  further  study  is  the  region  of  validity  of  the  plane- wave, 
integral  representation  for  the  electric  field  in  eq  (27).  Even  though  this  expression  gives  good 
results  when  analytically  continued  outside  a  spherical,  source-free  region  [39], [45],  the  validity 
of  this  analytical  continuation  has  not  been  rigorously  demonstrated.  This  question  has  practical 
consequences  for  field  uniformity,  stirring,  and  testing  since  rectangular  test  volumes  are  often 
used  in  practice. 
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APPENDIX  A.  RAY  THEORY 

The  mathematical  link  between  mode  theory  and  ray  theory  for  a  perfectly  conducting, 
rectangular  cavity  is  the  three-dimensional  Poisson  sum  formula  [19], [54].  This  formula  allows 
the  Green's  function  in  eq  (15)  to  be  converted  from  a  triple  sum  of  modes  to  a  triple  sum  of 
rays.  The  mathematical  details  are  fairly  involved  and  will  not  be  covered  here.  However,  the 
physical  interpretation  is  clearly  pictured  in  terms  of  multiple  images  as  shown  in  figure  19.  The 
source  is  a  z-directed  dipole,  and  the  multiple  images  represent  multiple  ray  bounces  in  they  =y' 
plane.  Similar  diagrams  could  be  generated  for  other  sources  and  locations. 

The  computation  of  the  field  at  a  point  in  the  cavity  is  tedious  because  of  the  triple  sum  of 
image  contributions.  In  fact  the  sum  is  not  convergent  for  some  frequencies  and  field  locations. 
This  has  to  be  the  case  because  the  mode  sum  in  eq  (15)  has  infinities  at  resonant  frequencies  of 
each  cavity  mode.  The  mode  representation  is  made  finite  for  imperfectly  conducting  walls  by 
introducing  a  finite  Q  (hence  the  resonant  frequencies  become  complex).  The  ray  sum  can  be 
made  finite  for  imperfectly  conducting  walls  by  introducing  a  reflection  coefficient  (which  has 
magnitude  less  than  one)  at  each  wall  bounce.  This  has  been  done  for  studying  the  field  buildup 
in  a  rectangular  cavity  when  the  source  is  a  turned-on  sinusoid  [55]. 

Multiple-image  theory  can  be  extended  to  include  the  effect  of  a  mechanical  stirrer.  Each 
image  cell  then  contains  an  image  of  the  mechanical  stirrer  with  appropriate  location  and 
orientation  as  shown  in  figure  20.  The  solution  of  the  large  boundary  value  problem  would  be 
extremely  difficult,  even  with  the  ray  tracing  approximation.  However,  the  multiple-image 
diagram  in  figure  20  can  be  used  to  provide  some  insight  into  stirrer  design.  The  goals  of  stirring 
are  to  randomize  the  field  and  to  eliminate  any  deterministic  component.  Another  way  to  state 
these  goals  is  to  minimize  the  ratio  of  unstirred  to  stirred  energy.  Unstirred  energy  arrives  at  the 
observation  without  interacting  with  the  stirrer.  An  example  is  (single-bounce)  ray  U  in  figure 
20.  A  stirring  strategy  then  is  to  design  the  stirrer  (or  stirrers)  to  eliminate  as  many  direct  rays  as 
possible.  The  conclusion  that  follows  is  that  the  stirrer(s)  must  be  comparable  to  chamber  size 
rather  than  just  comparable  to  a  wavelength.  This  conclusion  is  consistent  with  recent  chamber 
measurements  [8]. 


APPENDIX  B.  SHORT  ELECTRIC  DIPOLE 

Consider  a  short  electric  dipole  of  effective  length  L  oriented  in  the  z  direction  as  shown 
in  figure  21.  The  components,  Sra  and  Srp,  of  the  dipole  receiving  function  are  given  by 

S~=^p    and     S„=0,  (Bl) 

where  Rr  is  the  radiation  resistance.  In  eq  (Bl),  Sra  is  derived  by  dividing  the  induced  voltage  by 
twice  the  radiation  resistance  for  a  matched  load.  If  eq  (Bl)  is  substituted  into  eq  (63),  the 
angular  integration  can  be  carried  out  to  obtain 

E2L2 
<P>  =  ±^.  (B2) 

\2Rr 
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The  radiation  resistance  of  a  short  dipole  is  [20] 


Rr=^f.  (B3) 

r       3X2 

Substitution  of  eq  (B3)  into  eq  (B2)  yields  the  desired  final  result: 

1       T?2        32 

2  rj  Ak 

which  is  identical  to  eq  (66)  for  general  antennas.  The  polarization  mismatch  factor  of  one  half 
is  particularly  clear  for  the  electric  dipole  antenna  because  SrP  -  0 . 

The  idealized  electric-dipole  theory  in  this  appendix  does  not  apply  directly  to  electric- 
field  probes  because  most  probes  are  terminated  with  a  high  impedance  rather  than  a  matched 
load.  In  this  case  the  relevant  quantity  is  the  open-circuit  voltage  Vj  as  discussed  in  Section  4.2. 
Typically  only  the  field  magnitude  is  sensed: 

\K\*L\E.\  (B5) 

and  the  effective  length  is  determined  by  calibration  rather  than  theory.  Experimental 
comparisons  [8]  of  electric-field  components  with  field  probes  and  received  power  with  larger 
antennas  have  been  consistent  with  the  theory  in  Section  3 . 


APPENDIX  C.  SMALL  LOOP  ANTENNA 

The  other  electrically  small  antenna  of  practical  interest  is  the  small  loop.  For  a  small 
loop  of  area  A  centered  on  the  z  axis  in  the  xy  plane  as  shown  in  figure  22,  the  components  of  the 
receiving  function  are 


-  icofjA  sin  a 
~2~VRr 


Sra=0    and     Srp  =  -,UJr:       ■  (CO 


The  results  in  eq  (CI)  are  obtained  by:  (1)  determining  the  magnetic  flux  density  penetrating  the 
loop,  (2)  multiplying  by  -ico  to  determine  the  induced  voltage,  and  (3)  dividing  by  2Rr  to 
determine  the  current  induced  in  a  matched  load.  If  eq  (CI)  is  substituted  into  eq  (63),  the 
angular  integration  can  be  carried  out  to  obtain 

j-2      2      2  a  2 

\27j2Rr 
The  radiation  resistance  of  a  small  loop  is  [20] 
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_       2nr](kA\ 

A,    =  


5     \AJ 

Substitution  of  eq  (C3)  into  eq  (C2)  yields  the  desired  final  result: 


(C3) 


<Pr>  =  -^-—.  (C4) 

2  tj  An 

which  is  identical  to  eq  (66)  for  general  antennas  and  eq  (B4)  for  a  short  electric  dipole.  The 
polarization  mismatch  factor  of  V2  is  also  clear  for  a  small  loop  because  Sra  =  0 . 

APPENDIX  D.  CHAMBER  TIME  CONSTANT 

The  chamber  time  constant  is  of  interest  for  transient  applications  and  for  measurement 
of  Q  [7], [44].  Consider  first  the  case  where  a  cw  source  is  instantaneously  turned  off.  By 
equating  the  change  in  cavity  energy  Uto  the  negative  of  the  dissipated  power  times  a  time 
increment  dt,  we  obtain  the  differential  equation 

dU  =  -Pddt.  (Dl) 

Equation  (91)  can  be  used  to  replace  Pd  in  eq  (Dl): 

dU  =  -(coU/Q)dt  =  -  —  dt,  (D2) 

T 

where  the  time  constant  r  =  Q I  co  .  The  initial  condition  is  U=  Uq  at  /  =  0.  The  solution  of  eq 
(D2)  with  this  initial  condition  is 

C/  =  £/0exp(-f/r),     f>0.  (D3) 

The  closely  related  case  of  a  turned-on  (step  modulated)  cw  source  involves  the  same 
exponential  function  and  time  constant: 

£/  =  £/0[l-exp(-//r)],    />0.  (D4) 

The  average  energy  density  <W>  and  received  power  <Pr>  follow  the  same  exponential 
decay  or  buildup.  If  the  average  received  power  is  fit  to  an  exponential  curve,  the  inferred  time 
constant  rcan  be  used  to  measure  the  cavity  Q  [46], [56]: 

Q  =  cot.  (D5) 

This  method  has  the  advantages  that  it  does  not  require  a  knowledge  of  the  chamber  volume  and 
it  is  less  affected  by  antenna  efficiency  and  mismatch. 
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Table  1 .  Mode  properties  for  a  rectangular  cavity. 


Designation         Restrictions 


Nonvanishing  components  Number  of  modes 


TM/m„ 

/>l,/w>l,  n>\ 

t^x-,  t^y,  t^z,  tix,  tlz 

1  F'lmn 

t>\,m>\,  n>\ 

E>x,  &)>■>  tlx,  tiy,  tlz 

TM/m0 

l>\,  m>\,  n=0 

Ez,  HX,  Hy 

TEo/nn 

/=0,  m>\,  ri>\ 

EX,  Hy,  Hz 

TE/o„ 

l>\,  m=0,  ri>\ 

Ey,  Hx,  Hz 

Ni{k) 
N2(k) 
N3(k) 

N4(k) 
N5(k) 


LLJuLLLL 


Figure  1.  Cavity  of  volume  Fwith  perfectly  conducting  walls. 


39 


Figure  2.  Empty  rectangular  cavity  with  dimensions  axb  xc. 


Figure  3.  Source  current  distribution  J(r')  in  a  rectangular  cavity. 
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Figure  4.  Mode  number  (curve  1),  smooth  approximation  to  the  mode  number  (curve  2),  Weyl 
approximation  to  the  mode  number  (curve  3),  and  smooth  mode  density  (curve  4)  for 
the  NIST  reverberation  (2.74  m  x  3.05  m  x  4.57  m). 
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Figure  5.  Transmitting  antenna  in  a  reverberation  chamber  with  a  mechanical  stirrer. 
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►  y 


Figure  6.  Plane-wave  component  F(Q)  of  the  electric  field  with  wavenumber  k. 


Figure  7.  Unit  vectors,  s,  and  s2 ,  with  an  angular  separation  y. 


43 


Reverberation  Chamber 


Transmit 
Antenna 


Stirrer 


Microstrip 


Reference 
Antenna 


Figure  8.  Test  setup  for  immunity  measurement  of  a  microstrip  transmission  line  in  a 
reverberation  chamber. 
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Figure  9.  Comparison  of  theory  with  three  measurements  of  the  radiated  immunity  of  a 
microstrip  transmission  line  [40]  in  the  NIST  reverberation  chamber. 
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Figure  10.   Comparison  of  the  measured  probability  density  function  of  a  single  rectangular 
component  of  the  electric  field  [8]  with  theory  (Rayleigh  distribution). 
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Figure  11.   Comparison  of  the  measured  probability  density  function  of  the  total  electric  field  [8] 
with  theory  (chi  distribution  with  6  degrees  of  freedom). 
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Figure  12.   Comparison  of  the  measured  probability  density  function  of  received  power  [8]  with 
theory  (exponential  distribution). 
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Figure  13.  Plane- wave  reflection  from  an  imperfectly  conducting  wall  of  a  reverberation 
chamber. 
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Figure  14.   Comparison  of  Q  measured  by  power  ratio  (Qm:Loss)  and  decay  time  (Qm:TC)  with 
Q  calculated  from  eq  (97)  for  an  aluminum  cavity  [44].  The  theoretical  values  for 
wall  loss  (Ql)  and  receiving  antennas  (Q4)  are  also  shown. 
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Figure  15.   Comparison  of  Q  measured  by  power  ratio  (Qm:Loss)  and  decay  time  (Qm:TC)  with 
Q  calculated  from  eq  (97)  for  an  absorber-loaded  aluminum  cavity  [44].  The 
theoretical  values  for  wall  loss  (Ql),  absorption  by  salt-water  spheres  (Q2H),  and 
receiving  antennas  (Q4)  are  also  shown. 
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Figure  16.   Comparison  of  theory  with  three  measurements  of  the  radiated  emissions  of  a 
microstrip  transmission  line  [40]  in  the  NIST  reverberation  chamber. 
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Figure  17.  Equipment  under  test  (EUT)  radiating  fields  Et,  Ht  (emissions  measurement)  or 
illuminated  by  by  fields  Eu  Hi  (immunity  measurement). 


A 

p 

ZL 

1 

^^ 

1 

Figure  18.  Thevenin  equivalent  circuit  for  a  single  port  in  equipment  under  test. 
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Figure  19.  Multiple  images  for  a  z-directed  dipole  source  in  a  rectangular  cavity. 
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Figure  20.   Images  of  source  and  stirrer  in  a  rectangular  cavity.  The  (single-bounce)  ray  U  is  not 
affected  by  the  stirrer  (hence  contributes  to  unstirred  energy). 
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►  y 


Figure  21.  Dipole  antenna  illuminated  by  a  plane-wave  component  of  the  electric  field. 


►  y 


Figure  22.  Loop  antenna  illuminated  by  a  plane-wave  component  of  the  electric  field. 
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